Four-Dimensional Quantum Hall Effect in a Two-Dimensional Quasicrystal 
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One-dimensional (ID) quasicrystals exhibit physical phenomena associated with the 2D integer 
quantum Hall effect. Here we transcend dimensions, and show that a previously inaccessible phase of 
matter - the 4D integer quantum Hall effect - can be incorporated in a 2D quasicrystal. Correspond- 
ingly, our 2D model has a quantized charge-pump accommodated by an elaborate edge phenomena 
with protected level crossings. We propose experiments to observe these 4D phenomena, and gener- 
alize our results to a plethora of topologically equivalent quasicrystals. Thus, 2D quasicrystals may 
pave the way to the experimental study of 4D physics. 
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The three-dimensionality of our world has never pre- 
vented physicists from studying their theories in arbi- 
trary dimensions. In recent decades this has proven use- 
ful in mesoscopic physics, as many systems are well de- 
scribed by two, one and even zero dimensions. However, 
dimensions higher than three, which host many fascinat- 
ing phenomena, are seemingly out of reach. 

The uprising field of topological phases of matter is 
also dealing with systems of arbitrary dimension [HIHI- In 
this paradigm, each energy gap of a system is attributed 
an index, which is robust to continuous deformations. 
A nontrivial index is usually associated with interest- 
ing boundary phenomena, quantized response, and ex- 
otic quasiparticles. While nontrivial topological phases 
appear in any dimension [31 13] , the physical manifesta- 
tions are limited to ID [3] , 2D [3 |H] and 3D . 

An example for an intriguing topological phase, which 
is seemingly out of reach, is the 4D generalization of the 
2D integer quantum Hall effect (IQHE). In 2D, a uniform 
magnetic field creates Landau levels that are character- 
ized by the 1**' Chern number - the topological index that 
corresponds to the quantized Hall conductance [TUl HJ . 
In 4D, a uniform SU(2) Yang-Mills field results in gener- 
alized Landau levels [121 IE] • These levels are character- 
ized by the 2^^*^ Chern number - a topological index which 
corresponds to a quantized non-linear response |131 114j . 
Both the 2D and 4D IQHEs exhibit a variety of exotic 
strongly correlated phases when interactions are included 
[l3l I15j . Hence, the 4D IQHE, as well as lattice models 
with non- vanishing 2^^^^ Chern numbers, constantly at- 
tract theoretical attention [T^l [T7] . 

Recently, it was shown that ID quasicrystals (QCs) - 
non-periodic structures with long-range order - exhibit 
topological properties of the 2D IQHE [TH] . The bulk en- 
ergy spectrum of these ID QCs is gapped, and each gap 
is associated with a nontrivial I''' Chern number. This 
association relies on the fact that the long-range order 
harbors an additional degree-of-freedom in the form of a 
shift of the quasiperiodic order. Accordingly, boundary 
states traverse the gaps as a function of this shift. This 
property was observed in photonic QCs, and was utilized 
for an adiabatic pump of light }18j . Moreover, upon a 



deformation between QCs with different l*** Chern num- 
bers, a phase transition is expected to occur. Such a 
transition was also observed experimentally in photonic 
QCs [IH]. Generalizations to other ID symmetry classes, 
physical implementations, and QCs were discussed, see 
e.g. Refs. [M23]. 

In this work, we take a major step further, and present 
a 2D QC that exhibits topological properties of the 4D 
IQHE. Each gap in its energy spectrum is characterized 
by a nontrivial 2'^'^ Chern number, which implies: (i) 
quantized charge pumping with an underlying 4D sym- 
metry; (ii) gap-traversing edge states with protected level 
crossings; (iii) quantum phase transitions between topo- 
logically distinct QCs. Generalizations to other models 
and other types of QCs are discussed. We propose two 
experiments to measure the 2"^^ Chern number via charge 
pumping, and, thus, make 4D physics experimentally ac- 
cessible. 

We study a 2D tight-binding model of spin-^ particles 
that hop on a square lattice in the presence of a modu- 
lated on-site potential 
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Here Cx,y,a is the single-particle annihilation operator of 
a particle at site (x, y) with spin a; tx,ty are the hop- 
ping amplitudes in the x- and y-directions; and Ix, ly are 
the amplitudes of the on-site potentials, which are mod- 
ulated along X and y with modulation frequencies bx,hy, 
respectively [cf. Fig. [IJa)]. Lastly, the Hamiltonian de- 
pends on two shift parameters 4>x and (j)y We assume the 
modulation frequencies hx and by to be irrational, which 
makes the on-site modulations incommensurate with the 
lattice, and the model becomes quasiperiodic. 

The spectrum of H is gapped. Therefore, it may ex- 
hibit nontrivial topological indices. Note that H is nei- 
ther symmetric to time reversal nor to charge conjuga- 
tion. Therefore, according to conventional topological 
classifications [31 H] , the only apparent topological index 
that can be associated with it is the 1*^' Chern number. 
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FIG. 1. (a) An illustration of 7? [cf. Eq. a 2D square lat- 
tice with a modulated on-site potential that is spin dependent. 
The spin-up state, a = +, experiences the upper (blue) po- 
tential, and the spin-down state, a = —, experiences the lower 
(red) . For clarity, the potentials are vertically-displaced. The 
dots mark the underlying lattice sites, (b)-(d) A quantized 
charge-pumping along the a;-direction is achieved by scanning 
the shift parameters and (py in the presence of the modu- 
lation modifications By^ and -Buji/, respectively [cf. Eqs. Q 
and (|6|]. The effects on the potentials oi a — + (blue) and 
a — — (red) is illustrated (the thin (gray) line denotes the 
unmodified reference potential): (b) 4>x (or equivalently (j)y) 
shifts the cosine-potentials in opposite directions for opposite 
spins, (c) Byz makes the modulation frequency become spin 
dependent by + a Byz. (d) B^y shifts the cosine-potentials 
in opposite directions, but with an increasing shift along the 
j/-direction. 

However, for (jj^ — (l)y ^ 0, H is, m fact, time-reversal 
symmetric, and any H{(j>x, 4>y) is adiabatically connected 
to H{0,0). Hence, the I*'* Chern number of H vanishes, 
and H is seemingly trivial. 

Strikingly, H is attributed a nontrivial 2'^'* Chern num- 
ber, which, by definition, characterizes 4D systems. In 
order to obtain this result, let us consider H on a, toroidal 
geometry [21] . We introduce twisted boundary condi- 
tions along the x and y directions, parameterized by 
9x and 9y, respectively. For a given gap and given 
= i<Px,(l)v,0x,0y), we denote by P(0^) the projec- 
tion matrix on all the eigenstates of H{(j)^) with energies 
below this gap. We can now define, 

n(A\ ST '^o'PiS ( r, 9P d dP dP\ 



where eajs-ys is the antisymmetric tensor of rank-4. 

Formally, the 2'^'^ Chern number is defined by V = 
/ d'^(j)fj, C{(j)^) [25]. By its definition, V characterizes a 
4D family of Hamiltonians composed of all H{4>fj) with 
all possible values of 4>^, i.e. (jix and 4>y act as two addi- 
tional effective dimensions. The main result of this Let- 
ter is that even for a given (j)^, the gaps of _ff (0^) can be 
associated with a nontrivial integer 2"*^ Chern number, 

V = (2^)4C(0^)^O, (3) 

in the thermodynamic limit. Below, we justify Eq. ([s]) 
by showing that, for our model, C(0^) is essentially in- 
dependent of (j)^, and thus the integration over the four 
parameters is redundant. 

Beforehand, we present the physical implications of the 
nontrivial V. To do so, we apply to H the procedure of di- 
mensional extension that was introduced in Refs. [T51I21) . 
In this procedure, we interpret ipx and 0j, as momenta 
along two fictitious perpendicular coordinates w and 2, 
respectively. Now, the Hamiltonian H{(j)x, 4>y) is a single 
Fourier component of some ancestor 4D Hamiltonian T-L. 
By making the inverse Fourier transform, we obtain a 
Hamiltonian describing spin-^ particles hopping on a 4D 
hyper-cubic lattice 

n^Y. 4e''""'^^"^i/iCx+/i + h.c. , (4) 

where Cx = (cx,+ ,Cx,-) annihilates a spin-^ particle at 
site X = (x,?/, ji is summed over the unit vectors 

x,y,z and w, which connect nearest neighbors, and tjx = 
{tx, ty,Xx/2, Ay/2). The particles are coupled to a Yang- 
Mills gauge field a^(x) = [byyz + bxxw)^^. This vector 
potential describes a spin-polarized uniform SU(2) field. 
Such a field in known to generate a 4D IQHE with a 
nontrivial V [HI [13] . Notably, H is defined on a planar 
geometry in a Landau gauge, whereas previous analyses 
treated a spherical geometry in a symmetric gauge. 
Similar to the 2D IQHE, V has a physical manifestation 

in the form of a response function. Here the response is 
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quantized but non-linear: ja — V -^^^eap-isBij^Es [14j . 
where ja denotes the current density along the a- 
direction, $0 is the flux quantum. Eg is an electric field 
along the (5-dircction, and Bp~f is a magnetic field in the 
/37-plane. Note that B^^ couples to the electric charge, 
and affects equally both spin components. 

A direct observation of this response requires a 4D sys- 
tem. However, we can develop an analogue of Laugh- 
lin's pumping, which is manifested in the 2D QC. Let 

us consider the following two cases: jx — V -^^By^E^ 

2 

and jx — V j^B^yEz- Recall that the electric fields, E^ 
and Ez, can be generated by time-dependent Aharonov- 
Bohm fluxes, E^ S^5t*»(0 and E^ = ^dt<^z{t), 
where o is the lattice spacing, and and are the 
number of lattice sites along the w- and z-directions, re- 
spectively. Expressing Byz and B^y in H in Landau 
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gauge, and performing dimensional reduction, these fields 
enter the 2D model, H , through modified on-site terms, 

\x cos[2iT{ahxX + By^yu) + (l)x{t)] 

+ \yCos[2Ti{aby + By^)y + (j)y{t)], (5) 



where 4>x{i) 



2tt 



<^^{t) and (/)y(t) 
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Bwy = ByjyO? /(^Q and Byz = ByzO^/^o denote the cor- 
responding flux quanta per plaquette. Figures [ijb)-(d) 
illustrate the effects of these modifications. Now, fixing 
the chemical potential within a gap with a given V, an 
adiabatic scan of (l)^ or of 4>y from to 2n pumps charge 
along the x-direction, such that 



(t), and 



Qx ^^ByzNy , 

Qx = VeByjyNy , 



(6a) 
(6b) 



respectively, where Ny is the number of lattice sites along 
the 2/-direction. 

We can now propose experiments that measure V using 
Eq. For that purpose, take a 2D slab of the QC and 
connect metal leads to the edges of the a:-coordinate. Let 
us assume that the chemical potential of both the QC and 
the leads is placed in some gap of H . Then, one should 
measure the charge-flow during the scan of 4)x from to 
27r for different values of B 
we expect that V 



According to Eq. ( 6a 



-^dQx/dByz- Similarly, according 



to Eq. (6b), during the scan of 
charge flows and V = jj^dQx/dB^y 



while varying B^y, 
Remarkably, due 



to the 4D origin of our model, the measured V would be 
the same in both experiments. 

We have just seen that upon a scan of (j)x, charge may 
flow in the x-direction. Consequently, for an open geom- 
etry, in order to accommodate this charge transfer, edge 
states must appear and traverse the gaps as a function 
of (jjx- These states appear for infinitesimally small Byz, 
and hence appear also for Byz — 0. Figure [2] depicts the 
numerically obtained energy spectrum of if as a function 
of (l)x, for an open a:-coordinate, a periodic y-coordinate, 
tx^ty^ 1, Ix-^ly^ 1.8, x = v = f ..34, hx = (l + \/5)/2 
and by = 55/34, which approximates (1 -f a/5)/2 to or- 
der 10~^ |24| . The spectrum is invariant with respect 
to (f>y, and is depicted for (f)y = 0. As a function of 4>x, 
the spectrum has flat bands and gap-traversing bands. 
The flat bands correspond to bulk states, whereas the 
gap-traversing ones to edge state (see insets). The edge 
states are divided into four types: spin-up and spin-down 
states (blue and red), which are localized at either the left 
or right edge (opposite slopes). These edge states are a 
signature of the nontrivial V of our model. They can be 
measured in a way similar to the expreiments performed 
in ID photonic QCs [H]. 

Naively, one could postulate that opposite-spin modes 
that reside on the same edge could be gapped out by in- 
troducing spin-mixing terms. However, this is not the 




FIG. 2. The spectrum of // as a function of the shift param- 
eter (l>x [cf. Eq. ([1])], for an open x-coordinate and a periodic 
y-coordinate. The horizontal bands correspond to bulk states. 
The gap-traversing states with a = -\- (blue) or a = — (red) 
are edge states, which are locahzed at the left or right edges 
(see insets for typical wave functions). The crossings of cr = -|- 
and a = — bands are topologically protected. The values of 
the 2"^^ Chern number, V, associated with the large gaps are 
presented. 



case, and the edge modes and their crossings are topo- 
logically protected. In order to establish this protection, 
we decompose H into its spin and spatial constituents. 
The Hamiltonian iJ is a sum of two spin-components, 
where each spin-component is subject to two decoupled 
ID Harper models along the x- and y-directions. Both 
spin components experience the same modulation fre- 
quencies, hx and hy, but couple to the shift parame- 
ters, 4>x and 4>y, with an opposite sign. Therefore, each 
eigenstate of H{cj)x,(j)y) is a product of eigenstates of the 
Harper models in the x- and y-directions and a spin eigen- 
state. 

Recall that each gap of the Harper model is associated 
with a nontrivial I*'* Chern number, which corresponds to 
the number of boundary states that traverse the gap as a 
function of (j) dUl HH HH] • Accordingly, each band of the 
quasicrystalline Harper models in the x- and y-directions 
is associated with some spin-dependent Chern number, 
i^r^.a and i'ry,a, whcrc rx and Vy denote the correspond- 
ing bands, respectively [27] . In Fig. [2] the bands that tra- 
verse the gaps as a function of (px are composed of prod- 
ucts of bulk bands in the y-direction and boundary states 
in the a;-direction. Notably, due to the opposite coupling 
of the spin to (j)x and the Chern numbers of bands 
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of opposite spin have opposite signs, Vr^.+ = ~^r^,- and 
Vry,+ = ~'^ry,-- Therefore, the gaps are traversed by 
the same number of spin-up and spin-down bands, but 
with opposite slopes. Since opposite-spin bands are asso- 
ciated with opposite l*** Chern numbers, they cannot be 
gapped-out by spin-mixing terms, even if they cross at 
some value of (jjx- Otherwise, Vry,a would change contin- 
uously as a function of (px between i^ry,+ and i^ry.~- The 
level crossing is therefore protected. 

The described edge phenomena accounts for the 
charge-pumping appearing in Eq. When (px or (j)y 
are scanned, spin-up and spin-down states, which have 
opposite I''* Chern numbers, flow in opposite directions. 
In the absence of Byz and the two charge currents 
cancel. Applying By^ or B^jy causes a difference between 
the densities of spin-up and spin-down states, and thus a 
net charge is pumped ^4j. This corresponds to the idea 
of associating the 2""^ Chern number with pumping of l*** 
Chern numbers |28j . Remarkably, even for vanishing Byz 
and B^y, spin is pumped across the sample. Therefore, 
for recursive scans of (f)x, macroscopic spins accumulate 
at the boundaries. 

We turn, now, to establish Eq. ([s]). Let us evoke 
the definition of the f** Chern number of a band 
of the Harper model in the x-direction with spin cr, 

'^r^,o- = / d(l)xd0x Cr^.a{(t>x,dx), whcrC Cr^,a{<l>x, ^x) = 
^iT'r{Pr^A9<l>^Pr^,a,de^Pr^,a]), and Pr^,a{(i)x,Ox) is 

the projection matrix on the eigenstates of the Vxih. 
band |25j . A similar definition applies for ,tT- Let us 
denote by en{4>) the eigenenergies of the Harper model. 
The decomposition of H into x- and y-constituents makes 
its energy spectrum a Minkowski sum, -Bn^.n ,dr(0£c, (f'y) = 
^n:cio'4>x) + ^nyio-(t>v)- Accordingly, the states below each 
gap of En^^riy.a Can be decomposed into a sum over 
pairs of bands in the ID spectra, rx and Tj,, such that 
+ try < where fj, is the energy in the middle of the 
gap. Using the fact that the eigenfunctions of H are a 
product of the ID eigenfunctions, we obtain [22] 

V=$] l^r^^.l^ry^a^Q. (7) 

In a previous work [18], we have shown that, in the 
thermodynamic limit, C{(f>,9) becomes independent of 
and 9. Hence h'r,,^ = i'2'^)'^Cr^,a-{4'x,dx) and iyry,a = 
{2'KYCry^a{<t'y,&y)- This, Combined with Eq. Q, imme- 
diately implies Eg. ^ [29 ^ 

The invariance of C((j>, 9) of the Harper model with re- 
spect to (j) follows from the correspondence of to lattice 
translations. The modulated potential cos(27r6x -I- 0) im- 
plies that a translation by n lattice sites is equivalent to a 
shift of (j) by 27r(&rimod 1). Provided that h is irrational, 
the set of 0-shifts generated by all lattice translations is 
dense in the interval [0, 27r], in the thermodynamic limit. 
Consequently, any shift of cj) is equivalent, to arbitrary 
precision, to some lattice translation, which has no effect 



on C((/), 0). Turning back to i?, the invariance of C{(t)^) 
with respect to (j)x and (j)y results from the equivalence of 
their shifts to lattice translations in x and y, respectively, 
where the spin components are translated in opposite di- 
rections. 

Until now, our analysis used the spin-flip symmetry of 
H: a —> —a. Let us break this symmetry by making the 
modulation frequencies spin-dependent, i.e. bx — >■ bx,cr 
and by — )■ by a-- Now, a translation by n lattice sites, 
for example in the x-direction, is equivalent to spin- 
dependent shifts of (px by 27r((T62;^o-'T'mod 1). Ostensi- 
bly, (px is no longer equivalent to lattice translations. 
Nonetheless, if bx.+ and fex,- are mutually irrational, the 
Kronecker-Weyl theorem [30] implies that the space of 
spin-dependent shifts generated by lattice translations 
along X is dense over [0, 27r]^. In particular, the line 
created by shifts of (px is also densely covered. Conse- 
quently, here too, any shift of (j)x is equivalent, to arbi- 
trary precision, to some lattice translation. The same 
argumentation applies in the y-direction. Thus, Eq. ([S]) 
holds also in the absence of spin-flip symmetry. 

The above symmetry breaking may result from a U(l) 
gauge field in the 4D Hamiltonian H. In fact, any 
SU(2)xU(l) gauge transformation in 4D that respects 
Landau gauge keeps the system unchanged. After the 
dimensional reduction to 2D, such a transformation be- 
comes a general local transformation that may mix x, y, 
and spin, but keeps V unchanged. More generally, the 
properties of C(0p) implies that any unitary transforma- 
tion of the 2D Hamiltonian that does not depend on 0^ 
keeps C{(pfj_) independent of </)^. 

The predicted 2"^^ Chern number is not limited to a 2D 
QC composed of two Harper models. One can consider 
(i) placing the cosine modulations in the hopping terms, 
rather than in the on-site terms (off-diagonal Harper), 
and (ii) replacing each of the cosine modulations with 
a Fibonacci-like modulation. In ID, all these variants 
were shown to be topologically equivalent to the Harper 
model, namely they have the same distribution of 1*** 
Chern numbers [21 . Therefore, the gaps in such 2D QCs 
variants will have nontrivial V. 

To conclude, in this Letter we have presented a novel 
2D quasicrystalline model that is associated with the 
same topological index as the 4D IQHE - the 2"^^ 
Chern number. Correspondingly, our model exhibits 
an elaborate edge phenomena and a quantized charge- 
pump. We propose an experiment in which the charge is 
pumped through the system following a modification of 
the quasiperiodic modulation. Interestingly, while these 
modifications differ considerably, they lead to the same 
pumped charge. This equivalence may seem baffling from 
a 2D perspective, but follows from a simple symmetry 
of the corresponding non-linear response in 4D. Recent 
progress in controlling and engineering systems, such as 
optical lattices |31j and photonic crystals [ilJ? , makes our 
2D model seem experimentally feasible. Moreover, in- 
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teractions in such systems may lead to fractional qua- 
sicrystalline phases which are descendants of the exotic 
4D fractional quantum Hall effect [13]. Thus, our model 
serves as a porthole by which to access 4D physics. 
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SUPPLEMENTAL MATERIAL 



In the main text, we notice that the Hamiltonian H 
[cf. Eq. (1) of the main text] can be decomposed into 
Harper models in the x- and y-directions, for both spin 
species. Consequently, we infer that the 2"^^ Chern num- 
bers V of -ff can be expressed as multiplications of the 
l*** Chern numbers of the constituting Harper models 
[cf. Eq. (7) of the main text]. In this Supplemental Ma- 
terial we prove this preposition. 

Recall that for a given gap in the spectrum of 
we define the projection matrix P{4>^) on all the eigen- 
states of H{(t)^) with energies below this gap. Similar 
to iJ, P = Scr=± ^o-' where projects on the (T-spin 
states, and thus P+P- = 0. Additionally, we note that 
the states below the gap of the 2D model, H , are com- 
posed of pairs of ID bands that satisfy e^^ + er„ < 
where Er is an eigenenergy of the rth band of the Harper 
model, and /i is an energy in the middle of the 2D gap. 
Therefore, 



where is the antisymmetric tensor of rank-2, and 
C{(j),9) is the 1*** Chern form. The 2"^^ Chern number 
corresponding to C is 



(1.3) 



d(j)yd9y Cr^.a = J^i 



where z^rj^o- and i'ry,a are the I''* Chern numbers associ- 
ated with Pr^,a{4>x,0x) and Pry,a{<t)y,Oy), respectively. 
V is the 2""^ Chern number of a single term in the sum 



(I.ll, i.e. with given r^, Vy and a. Summing over all 
terms, we obtain the total 2'^'^ Chern number, V, given 
by Eq. (7) of the main text. 



n^p,) = E E 



(1.1) 



where Pr^.a{4>x^ 9x) is the projection matrix on the eigen- 
states of the ra;th band, and similarly for Pry.a{(t>y,9y). 
Note that the elements of this sum are orthogonal pro- 
jectors. 

In the main text, we introduced the 2"*^ Chern form, 
C(0^) [cf. Eq. (2) of the main text], and the 2"'^ Chern 
number, V — J C(0^). Since V is a Z- index, a 2"^^ 
Chern number can also be assigned to each element of 
the sum in Eq. ( |I.1[ ). It suffices, therefore, to focus on a 
single such element, and show that its 2"*^ Chern number 
is a product of its corresponding 1"* Chern numbers. 
Let us denote such an element by P{4>^i) = 

and its corresponding 



P, 



Pr^,a{4>x,Qx, 

2"*^ Chern form by C{4)^). For brevity, we de- 
note dP/d(j)a by p["l, which makes C{(j)^) = 
e„^^5Tr(PP["lp['51ppWp[''l)/(-87r2). According to the 
decomposition of P into x- and y-parts, C involves two 
types of terms: 



(\\ Tr(P P'"J P P'^' P 61) P P"^' P P'' 



(ii) Tr(P,^,,P^Sp[SP, 



p p[7] p['5] \ 



The type-(i) terms vanish, since for any projector, 
P, the following applies: PP^'^^P = 0, where p["l 
is a derivative of P with respect to some variable 
a. Recall that for any two matrices, Ti{A (S) B) ~ 
Tr(A)Tr(P). This makes the type-(ii) terms be- 
come Tr(P,^,,PJSpJS)Tr(P,^,,Pi;;'!<,pJ^l,<,). Substitut- 
ing this into C, we obtain 



C(0m) = gTr(P.,,.pH,p[S) • gTr(P.^,.pW,pW J 

= Cr^,ai4'x,0x) Cry,cr{4>y,(^y) , (1-2) 



